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Abstract 

Pulse-like  acoustic  signals  aFe, transmitted  from  an  acoustic  source  near  Oahu  to  seven 
receivers  off  the  west  coast  of  the  United  States  for  a  124-day  period  in  1988.  Acoustic 
travel-time  oscillations'iare  observed  in  the  received  signal  at  periods  between  15  and  23 
hours,  which  )*fe  caused  by  barotropic  (or  first  or  second  mode  baroclinic)  fluctuations 
in  the  ocean.  It  is  shown  that  these  fluctuations  cannot  be  local  processes  isolated  to 
either  the  source  or  to  the  receivers.  It  is  further  shown  that  resonant  barotropic  gravity 
wave  modes  (Platzman  et  al.,  1981)  are  not  consistent  with  the  data.  The  cause  of  these 
fluctuations  remains  unresolved,  but  the  data  and  other  oceanographic  measurements  put 
many  constraints  on  the  process  causing  these  fluctuations. 


Thesis  Supervisor:  Dr.  John  L.  Spiesberger 
Woods  Hole  Oceanographic  Institution 
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Chapter  1 


Introduction 

1.1  Motivation 

Acoustic  tomography  has  been  successfully  used  to  analyze  mesoscale  (order  lOO-lOOOkm) 
features  of  the  ocean  (Spiesberger  et  al.,  1983;  DeFerrari  and  Nguyen,  1986;  and  others). 
For  these  experiments,  other  synoptic  methods  of  extracting  information  about  mesoscale 
features  are  available  (i.e.  XBT’s,  CTD’s,  and  current  meters).  Yet  when  considering 
synoptic  measurements  on  a  basin-scale,  the  cost  of  deploying  such  equipment  becomes 
prohibitive.  In  the  past,  analytic  and  numerical  methods  have  been  relied  upon  to  provide 
dynamic  models  of  these  features.  An  acoustic  tomography  array  can,  at  reasonable 
cost,  be  used  to  obtain  synoptic  measurements  of  barotropic,  basin-scale  oceanographic 
phenomena  (Spiesberger  et  al.,  1989b). 

A  four  month  basin-scale  tomography  experiment  was  performed  in  1988  in  the  north¬ 
east  Pacific  where  a  bottom-mounted  source  near  Oahu  transmitted  a  continuous  (133Hz) 
acoustic  signal  to  seven  bottom-mounted  receivers  off  the  west  coast  of  the  United  States. 
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Signal  processing  converted  the  signal  at  the  receiver  into  an  equivalent  signal  which  would 
be  received  if  the  source  transmitted  a  pulse-like  signal  of  GO  ms  duration  at  184.2  s  inter¬ 
vals.  Changes  in  acoustic  travel-time  are  measured  by  estimating  the  phase  change  of  the 
received  signal  of  adjacent  records.  This  experiment  follows  up  two  earlier  experiments 
(in  1983  and  1987)  in  which  a  source  near  Oahu  transmitted  to  one  receiver  olT  the  coast 
of  northern  California  for  5  and  21  days,  respectively. 

These  earlier  experiments  yield  the  following  results.  Travel-time  oscillations  at  tidal 
frequencies  are  observed,  and  are  caused  by  barotropic  tidal  currents.  In  addition,  sig¬ 
nificant  oscillations  are  observed  at  many  non-tidal  frequencies.  Specifically,  at  periods 
between  15  and  23  hours,  oscillations  are  observed  having  rms  variations  of  about  10  ms 
which  vary  in  amplitude  and  frequency  on  about  a  10-day  time  scale,  indicating  that  they 
may  be  intermittently  excited.  Since  data  were  only  available  at  one  receiver,  it  has  not 
been  possible  to  match  these  oscillations  to  an  oceanographic  process  (Spiesberger  et  al., 
1989b) 

With  results  from  the  1988  experiment,  presented  in  this  thesis,  we  examine  four 
hypotheses  in  a  search  for  the  oceanographic  source  of  these  oscillations.  The  hypotheses 
are: 

•  The  ocean  fluctuates  only  near  one  receiver. 

•  The  ocean  fluctuates  independently  at  each  receiver. 

•  The  ocean  fluctuates  only  near  the  source. 

•  The  ocean  fluctuates  at  basin-scales  as  a  result  of  resonant  barotropic  gravity  wave 
modes. 
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We  show  each  of  these  hypotheses  to  he  unlikely,  and  are  not  able  to  provide  a  theory 
for  the  cause  of  the  observed  travel-time  oscillations  which  is  consistent  with  both  the 
tomographic  data  and  other  measurements  of  the  ocean. 

1.2  Overivew 

The  thesis  is  organized  as  follows  Chapter  2  contains  a  brief  discussion  of  acoustic 
propagation  in  the  ocean  as  it  relates  to  acoustic  tomography.  Chapter  3  discusses  the 
experiment  and  signal  processing.  Chapter  4  examines  each  of  the  four  hypotheses  pro¬ 
posed  to  explain  the  cause  of  the  acoustic  travel-time  oscillations.  Chapter  5  concludes 
with  a  discussion  of  the  constraints  placed  on  the  problem  by  both  the  tomographic  data 
and  other  measurements,  and  suggests  other  avenues  of  research  for  this  problem. 
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Chapter  2 


Acoustic  Propagation 


Sound  travels  at  about  1500  m-s-1  in  water,  which  is  about  four  and  one-half  times  that 
of  air.  Sound  speed  increases  with  increasing  temperature,  increasing  pressure,  or,  to  a 
much  smaller  extent,  increasing  salinity.  At  middle  and  equatorial  latitudes,  the  speed  of 
sound  initially  decreases  with  depth  because  of  cooling,  and  then  increases  in  the  deep 
ocean  due  to  great  pressure,  resulting  in  a  sound  speed  minimum  at  about  1  km  depth 
(Figure  2-1).  Application  of  Snell’s  law  of  refraction  shows  that  sound  is  refracted  towards 
regions  of  lowest  sound  speed  (Officer,  1958).  The  ocean  thus  has  a  sound  channel  (Ewing 
and  Worzel,  1948)  which  focuses  acoustic  energy  away  from  the  surface  and  bottom  where 
significant  attenuation  and  scattering  of  the  signal  would  otherwise  occur  (Figure  2-1). 
This  fact,  coupled  with  remarkable  water  transparency  to  sound  at  frequencies  below 
about  300  Hz  (Urick,  1983),  allows  detection  of  sound  propagated  over  thousands  of 
kilometers. 

At  frequencies  above  about  100  Hz,  propagation  of  sound  can  be  modeled  with  ray 
theory,  an  approximate  solution  of  the  acoustic  wave  equation,  since  the  acoustic  wave- 
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Figure  2-1:  Vertical  profile  of  sound  speed  in  the  northeast  Pacific  half-way  between  Oahu 
and  California  (left).  Twelve  acoustic  ra.y  paths  originating  from  a  source  on  the  sound 
channel  axis  for  a  200  km  section  of  ocean  (right). 


length  (about  15  m  at  100  Hz)  is  much  less  than  the  vertical  scale  of  the  channel  (about  1 
km).  Ray  theory  predicts  that  an  acoustic  pulse  emitted  from  a  submerged  source  travels 
to  a  distant  receiver  along  ray  paths  (multipaths)  which  have  different  travel-times  in  ac¬ 
cordance  with  Snell’s  law  of  refraction  (Officer,  1958).  Each  ray  path  samples  a  different 
vertical  slice  of  the  ocean,  based  on  the  departure  angle  (from  the  horizontal)  of  the  ray 
from  the  source.  Thus  the  travel-times  of  different  multipaths  can  be  used  in  an  inverse 
problem  to  solve  for  baroclinic  fluctuations  in  the  ocean.  Alternatively,  travel-times  from 
all  multipaths  can  be  combined  in  an  inverse  problem  to  solve  for  barotropic  fluctuations 
in  the  ocean. 


To  first  order,  fluctuations  in  acoustic  travel-time  can  be  computed  by  assuming  that 
the  ray  path  does  not  move  in  the  presence  of  fluctuations  (Munk  and  Wunsch,  1979). 
Perturbations  of  acoustic  travel-time,  6rt,  for  the  ith  ray  path,  f,,  are  primarily  related 
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to  in-sit u  fluctuations  of  temperature,  SO,  hydrostatic  pressure,  pyy.  and  current.  ?7,  ac¬ 
cording  to, 


S-r.it) 


S0(s,1) 

c'i(sJ) 


ds 


V(  f,  t ) 


ds 


u( s,t)  ■  ds 
<"o('S0 


(2.1) 


where  the  ith  ray  path  has  coordinates  parameterized  by  .7  (with  differential  elements  ds 
and  ds),  and  t  is  time.  The  reference  sound  speed  is  ca(s,t ),  and  that  portion  of  the 
sea-surface  displacement  related  to  hydrostatic  fluctuations  is  y.  The  water  density  is  p 
and  the  gravitational  acceleration  is  y.  The  constants  are, 

a  %  4.72  x  IQ-3  kin(°C  -s)_l 


/j  ss  1.7  X  10-2  s_1, 

(Spiesberger  et  al.,  1989b). 

Suppose  a  ray  path  travels  a  length  L  (km)  in  the  presence  of  a  uniform  oceanic  fluc¬ 
tuation  where  current  (parallel  to  the  ray  path)  is  u  (tn/s),  or  where  the  temperature 
fluctuation  is  SO  (°C),  or  where  hydrostatic  pressure  changes  are  caused  by  sea-surface 
displacement  y  (in).  Then  if  currents,  temperature,  or  sea-surface  oscillations  are  indi¬ 
vidually  responsible,  the  travel-time  oscillation  (s)  is  approximately, 


St  -4.44  x  10  4  L  u 

(2.2) 

St  %  -2.1  x  10"3  L  SO 

(2.3) 

St  rs  -7.55  x  10-6  L  y, 

(2.4) 

(Spiesberger  et  al.,  1989b). 


13 


Chapter  3 


The  Experiment  and  Signal 
Processing 

3.1  Experimental  Layout 

A  continuous,  phase-modulated  signal  (center  frequency  133  Hz.  bandwidth  17  Hz)  is 
transmitted  between  6  January  and  9  May  1988  (121  days)  from  a.  bottom-mounted  source 
off  the  north  coast  of  Oahu  (183  m  depth)  to  seven  bottom-mounted  receivers  ofT  the  west 
coast  of  the  United  States  (Figure  3-1).  The  receivers  are  about  3000  to  4000  km  from 
the  source  and  have  acoustic  travel-times  of  about  30  to  40  min  respectively.  The  source 
is  within  30  m  of  the  sources  used  for  the  1983  and  1987  experiments  (Spiesberger  et  al., 
1989b). 
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Figure  3-1:  The  1988  acoustic  tomography  experiment  consists  of  a  bottom-mounted 
source  (S)  located  near  Kaneohe  Bay,  Oahu  and  seven  bottom-mounted  receivers  (Rl  to 
R7)  whose  approximate  positions  are  shown.  The  distances  between  the  source  and  the 
receivers  are  about  3000  to  4000  km. 


3.2  Signal  Processing 


[n  order  to  raise  signal- to- noise  ratio  above  a  detectable  threshold,  the  transmitted  signal 
is  phase-modulated  every  eight  cycles  using  a  511-digit  pseudo  random  code.  A  matched 
filter  at  each  receiver  compresses  the  energy  of  the  30.7  s  cycle  into  that  of  an  equivalent 
transmitted  pulse  of  60  ms  duration.  Six  such  cycles  are  averaged  together,  yielding  travel- 
time  records  at  18-1.2  s  intervals,  further  increasing  signal-to-noise  ratio.  This  processing 
is  the  same  as  that  used  in  prior  experiments  (Spiesberger  et  al.,  1989a;  Metzger,  1983). 

At  each  receiver,  we  measure  changes  in  acoustic  phase,  60,  at  intervals  of  184.2  s 
using  a  maximum  likelihood  estimator  which  computes  the  phase  of  the  complex  cross¬ 
correlation  of  adjacent  travel-time  records  (Spiesberger  et  al.,  1989a).  This  estimator 
is  optimally  designed  to  detect  barotropic  fluctuations  as  well  as  first  or  second  mode 
baroclinic  fluctuations,  and  is  designed  to  attenuate  signals  from  higher  vertical  modes. 

If  the  acoustic  travel-time  fluctuates  at  period  T,  then  the  change  in  travel-time  is 
related  to  the  change  in  acoustic  phase  by. 


6t  - 


T6(j> 
u>cAt ' 


(3.1) 


where  the  source  frequency,  uc,  equals  2m - 133  rad-s-1  and  the  time  between  travel-time 
records,  At,  equals  184.2  s  (Spiesberger  et  al.,  1989a).  Travel-time  changes  are  related  to 
oceanographic  fluctuations  according  to  Equation  2.1. 

We  are  interested  in  oceanographic  fluctuations  other  than  internal  waves  of  higher 
mode  number  (>  2).  The  acoustic  signals  associated  with  these  internal  waves  are  con¬ 
sidered  to  be  a  portion  of  the  “noise.”  The  precision  of  the  travel-time  estimates  ( 135/xs 
at  At  =  184.2  s)  is  limited  by  the  higher  mode  internal  waves  rather  than  by  clock  errors 
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at  the  source  and  receiver  or  by  acoustic  noise  (Spiesherger  et  al.,  19S!)a). 
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Chapter  4 


Hypotheses 


We  consider  four  hypotheses  to  explain  measured  fluctuation  in  acoustic  travel-time  at 
periods  between  the  semi-diurnal  and  diurnal  tides.  They  arc: 

1.  The  o'-ean  fluctuates  only  near  one  receiver  (receiver  R3  of  Figure  3-1).  Previous 
papers  only  analyzed  data  from  R3  (Spiesberger  et,  al.,  1989a, b).  In  section  4.1,  we 
find  this  hypothesis  to  be  false. 

2.  The  ocean  undergoes  independent  fluctuations  near  each  receiver.  In  section  4.2,  we 
discuss  evidence  which  weighs  against  this  hypothesis. 

3.  The  ocean  fluctuates  only  near  the  acoustic  source.  In  section  4.3,  we  discuss  evi¬ 
dence  which  disagrees  with  this  hypothesis. 

4.  The  ocean  fluctuates  at  basin-scales.  In  section  4.4,  we  present  evidence  that  the 
oscillations  are  not  resonant  barotropic  gravity  wave  modes  of  the  world’s  oceans. 
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4.1  Hypothesis  1:  Ocean  Fluctuations  Localized  to  R3 


All  sovon  receivers  exhibit  prominent  oscillations  in  acoustic  travel-time  between  the  semi¬ 
diurnal  and  diurnal  tidal  periods  (Figures  1-1  and  1-2).  Therefore  this  hypothesis  is  not 
true. 


4.2  Hypothesis  2:  Independent  Ocean  Fluctuations  Near 
Each  Receiver 

We  first  check  if  acoustic  travel-time  oscillations  of  the  same  frequency  are  simultaneously 
observed  at  all  receivers.  The  121-day  experiment  is  divided  into  twelve  consecutive 
10-day  segments.  Periodograms  are  computed  for  each  receiver  which  collected  data 
during  the  segment.  We  find  that  the  oscillations  between  periods  of  15  and  23  hours  are 
simultaneously  observed  by  most  receivers  (Figure  4-3).  For  example,  between  yeardays 
70  and  79,  all  seven  receivers  exhibit  an  oscillation  at  periods  between  16.1  and  17.5 
hours  and  exhibit  an  oscillation  at  periods  between  19.0  and  20.9  hours  (Figure  4-3).  We 
conclude  that  most  receivers  simultaneously  observe  travel-time  oscillations  at  the  same 
frequency. 

We  believe  the  above  data  strongly  weighs  against  hypothesis  2.  Further  evidence 
against  this  hypothesis  is  presented  next. 

The  rms  variations  of  acoustic  travel-time  are  about  10  ms  for  each  peak  observed 
in  the  10-day  periodograms  at  R1  to  R7  (Figure  4-3).  In  1983  and  1987,  acoustic  os- 
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Figure  4-1:  Power  spectral  densities  of  a  42-day  time  series  (yeardays  41  to  83,  1988) 
of  phase  differences  at  each  receiver.  A  Bartlett  window  is  applied  to  the  data.  The 
mean  spectral  level  due  to  random  walk  noise  (Spiesberger  et  al.,  1989a)  is  at  or  below 
1  x  10-5cycles2/cpd.  The  vertical  solid  lines  mark  prominent  tidal  periods.  The  vertical 
dashed  lines  mark  periods  of  all  resonant  barotropic  gravity  modes  with  periods  between 
15  and  23  hours  (Platzman  et  al.,  1981). 
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Figure  4-2:  RMS  travel-times  (ms)  are  shown  for  the  peaks  in  Figure  4-1  with  periods 
between  15  and  23  hours  (vertical  axis).  The  arrows  (right  axis)  mark  the  periods  of  all 
resonant  barotropic  gravity  modes  with  periods  between  15  and  23  hours  (Platzman  et 
al.,  1081). 
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Figure  4-3:  The  number  of  receivers  exhibiting  spectral  peaks  within  a  given  range  of 
periods  is  shown  for  each  10-day  segment  of  the  experiment.  The  number  of  receivers  for 
which  data  arc  available  for  each  segment  is  noted  in  parenthesis  in  the  left  column.  In 
each  case,  the  spectra  are  computed  using  a  Bartlett  data  window.  The  periods  listed 
(top)  are  the  eigenfrcquencies  of  the  periodogram.  The  spacing  between  eigenfrequencies 
is  approximately  half  the  resolution  bandwidth  of  the  periodogram  (Harris,  1978). 


filiations  near  a  15  hour  period  had  similar  size  (Spiesberger  et  al.,  1989a,b).  Applying 
1'quations  2.2  and  2.3.  we  calculate  how  large  a  barotropic  current  or  how  large  a  temper¬ 
ature  perturbation  must  be  to  yield  an  rms  variation  of  10  ms.  If  the  ocean  fluctuations 
are  confined  within  a  50  km  region  of  each  receiver,  a  uniform  current  of  0.6  m  s-1  or  a 
uniform  temperature  perturbation  of  0.132°C  is  required.  Magnitudes  of  the  perturbation 
in  current  and  temperature  are  inversely  proportional  to  the  horizontal  scale  of  the  fluctu¬ 
ation.  If  the  fluctuation  is  confined  to  25  km,  then  the  associated  current  or  temperature 
perturbation  would  be  twice  the  above  valTies.  The  magnitudes  of  these  perturbations  arc 
implausible. 

We  therefore  find  it  unlikely  that  independent  ocean  fluctuations  near  each  receiver 
give  rise  to  the  observed  oscillations  in  acoustic  travel-time  at  periods  between  the  semi¬ 
diurnal  and  diurnal  tides. 


4.3  Hypothesis  3:  Oceanic  Fluctuation  Localized  Near  the 
Acoustic  Source 

We  present  two  results  which  indicate  that  the  acoustic  oscillations  with  periods  between 
15  and  23  hours  are  not  caused  by  ocean  fluctuations  confined  to  the  region  of  the  acoustic 
source.  These  results  are  presented  separately  in  subsections  4.3.1  and  4.3.2  below. 

4.3.1  Island  Trapped  Waves 

A  model  for  the  trapping  of  super-inertial  internal  waves  around  islands  was  proposed 
by  Wunsch  (1972)  to  explain  thermal  fluctuations  near  Bermuda.  In  this  model  the 
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ocean  (constant  depth,  rigid  lid),  containing  a  cylindrical  island,  is  realistically  stratified 
allowing  for  free  internal  waves.  Internal  wave  scattering  exhibits  pseudo-resonances  if  the 
period  of  the  incident  wave  matches  the  free  wave  travel-time  around  the  cylinder  (island). 
When  this  model  is  applied  to  Oahu  (Appendix  A),  no  significant  pseudo-resonances  are 
observed  for  wave  periods  between  10  and  25  hours. 

A  second  mechanism  for  the  generation  of  island  trapped  waves,  refraction  (resulting 
in  "trapped-leaky”  waves),  is  not  likely  to  occur  around  Oahu  due  to  the  narrow  shelf 
width  and  complicated  geometry  of  the  island  (Luther,  1985:  Longuet-IIiggins,  1967; 
Summerfield,  1972).  Thus,  available  theoretical  models  fail  to  provide  a  mechanism  for 
the  generation  of  island  trapped  waves  around  Oahu. 

Even  though  we  find  no  t  heoretical  justification  for  island  trapped  waves,  these  theories 
are  only  approximations  which  contain  some  cpiestionable  assumptions,  and  such  waves 
may  exist.  For  example,  Luther  (1985)  has  analyzed  10-  and  14-year  records  of  sea-level 
at  two  stations  on  opposite  sides  of  Oahu  (Mokuoloe  and  Honolulu).  His  analysis  reveals 
peaks  in  the  sea  level  spectra  from  each  station  at  periods  of  17.0  and  20.0  hours  with 
amplitudes  of  about  0.5  cm.  Cross-spectral  phase  relationships  suggest  that  these  spectral 
peaks  could  result  from  island  trapped  waves  with  horizontal  wavelengths  equal  to  1  or  | 
the  circumference  of  Oahu.  Patzert  and  Wyrtki  ( 1974)  have  observed  clockwise  currents  of 
about  5  cm/s  around  Oahu,  providing  further  support  for  the  existence  of  island  trapped 
waves. 

If  there  is  an  island  trapped  wave  of  the  type  hypothesized  by  Luther,  this  is  not 
the  wave  observed  in  the  tomographic  measurements.  Consider  the  transmission  paths 
to  R3  and  R4  which  leave  the  source  at  bearing  angles  differing  by  about  10°  (Figure 
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1-1).  If  a  15.5  hour  period  island  trapped  wave  has  a  horizontal  wavelength  of  i  the 
island  circumference,  then  the  predicted  phase  lag  between  acoustic  signals  at  R3  and 
R4  is  5°  (Figure  1-1).  Cross-spectral  phases  at  15.5  hour  period  (at  R3  and  R  1 )  have 
values  of  157°  ±  15°  (Figure  4-5  and  Appendix  B).  Calculations  of  cross-spectral  phases 
at  other  receiver  combinations  and  periods  yield  results  which  are  also  inconsistent  with 
the  island  trapped  wave  hypothesis  of  low  horizontal  mode  number  (Figure  4-5).  If  the 
horizontal  mode  number  of  the  island  trapped  wave  is  64,  then  the  predicted  phase  lag  at 
R3  and  R  1  would  match  the  observed  phase  lag.  We  do  not  believe  t  hat  realistic  coastlines 
and  bathymetric  features  could  support  such  an  island  trapped  wave.  'Flic  tomographic 
measurements  are  therefore  inconsistent  with  the  island  trapped  wave  hypothesis. 


4.3.2  Other  Near-Source  Fluctuations 

Unrealistically  large  fluctuations  of  current  and  temperature  are  required  to  yield  acoustic 
travel-time  oscillations  of  10  ms  rms  (Section  4.2).  We  conclude  that  the  tomographic 
measurements  are  inconsistent  with  hypothesis  3. 

4.4  Hypothesis  4:  Resonant  Barotropic  Gravity  Wave  Modes 

We  investigate  the  hypothesis  that  the  acoustic  oscillations  are  caused  by  resonant  barotropic 
gravity  wave  modes. 
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to  R3 


Figure  4-4:  The  relationship  is  shown  between  the  phase  lag,  n'F,  of  an  island  trapped 
wave  (horizontal  mode  n)  at  the  wave’s  intersection  with  two  ray  paths  and  the  difference 
in  bearing  angle  of  the  ray  paths,  0.  The  tomographic  signal  generated  by  an  island 
trapped  wave  will  have  a  phase  lag  between  the  two  receivers  equal  to  about  n<f>,  where 
n<P  =  n^Q.  The  distances  ri  and  r 2  are  about  40  and  10  km  respectively,  and  0  is  about 
10°.  The  angles  in  this  figure  are  enlarged  for  clarity. 
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RECEIVER 


Figure  4-5:  Phase  lags  (degrees)  are  shown  for  periods  at  which  coherency  peaks  occur 
in  the  cross-spectral  estimates  of  42-day  time  series  of  phase  differences  at  each  receiver 
referenced  to  R3.  Each  time  series  has  the  same  reference  time  and  runs  from  yearday  41 
to  83,  1988.  The  error  bars  show  the  range  of  periods  over  which  the  phase  shift  remains 
within  ±15°  of  its  value  at  the  coherency  peak  (Appendix  13).  The  arrows  (right  axis) 
mark  the  periods  of  all  resonant  barotropic  gravity  modes  with  periods  between  15  and 
23  hours  (Platzman  et  al.,  1981). 


4.4.1  Resonant  Shallow  Water  Gravity  Modes 


Platzman  (1978)  and  Platzman  et  al.  ( 1981 )  developed  a  model  identifying  normal  modes 
of  the  world  ocean.  In  this  model,  the  primitive  equations  are  applied  to  a  finite  element 
discretization  of  the  world  ocean  (about  900  km  between  nodes)  yielding  eigenfrequencies 
of  resonant  barotropic  vorticity  and  gravity  modes.  Corresponding  to  each  eigenfrequencv 
is  a  set  of  eigenvalues  describing  the  sea.  level  and  velocity  structure  of  each  mode.  Platz- 
man  et  al.  (1981)  find  eight  gravity  modes  at  periods  between  15  and  23  hours  (Figures 
1-1,  1-2,  and  4-5),  two  of  which  (15.5  and  21.2  hours)  are  dominant  in  the  Pacific  Ocean. 

In  our  study,  we  use  the  same  programs  developed  by  Platzman  to  identify  resonant, 
barotropic  gravity  modes  for  the  Pacific  Ocean  only  (Miller  et  al.,  1988).  Using  the 
same  grid  as  Platzman  (with  minor  variations  at  boundary  nodes),  we  find  three  resonant 
modes  at  periods  between  15  and  23  hours  (15.06,  18.46,  and  21.66  hours).  Considering 
the  altered  boundary  nodes,  this  is  consistent  with  the  Platzman ’s  results. 

For  each  mode,  we  compute  the  predicted  acoustic  travel-time  oscillation,  at  receiver 
i.  using  Equat  ion  2.1  and  obtain  results  of  the  form, 

fiTt(t)  =  aa,  sin(±a t  +  <&  +  (),  (4.1) 

where  a,  and  </>,  are  the  amplitude  and  phase  of  the  acoustic  travel-time  computed  from 
the  model,  a  is  frequency  and  t  is  time.  Since  the  model  magnitude  and  phase  of  the  sea- 
surface  and  velocity  are  based  on  arbitrary  references,  scaling  factors  a  and  are  required 
when  comparing  the  predicted  and  observed  acoustic  travel-time  oscillations.  For  each 
mode,  C  is  chosen  such  that  the  phase  of  the  predicted  acoustic  travel-time  at  R3  (the 
reference  receiver)  is  0°,  and  o  is  chosen  such  that  for  the  amplitude  of  the  predicted 
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acoustic  travel-time  matches  that  of  R3.  We  estimate  the  accuracy  of  our  model  travel 
times  to  be  about  20%  in  magnitude  and  ±15°  in  phase. 

A  comparison  of  the  observed  and  predicted  acoustic  travel-times  at  each  receiver  are 
shown  in  Tables  4.1,  4.2,  and  4.3  for  each  of  the  three  modes.  The  model  phase  at  each 
mode  deviates  significantly  from  the  observed  phase.  For  example,  at  the  15.06  hour 
mode  there  is  a  consistent  —160°  discrepancy  between  the  model  and  the  observations. 
We  conclude  that  the  observed  travel-time  oscillations  are  not  caused  by  a  single  mode. 

We  investigate  the  case  where  two  modes  occur  simultaneously  at  each  modal  fre¬ 
quency  (clockwise  and  counterclockwise).  We  find  no  combination  of  these  modes  at  any 
of  the  three  modal  frequencies  yield  the  observed  phase  and  magnitude  of  acoustic  travel- 
time  oscillations  at  each  receiver  (Appendix  C).  Wc  therefore  conclude  that  the  two-mode 
hypothesis  is  not  consistent  with  the  observations. 

4.4.2  Comparison  With  Other  Measurements 

For  each  resonant  barotropic  shallow  water  gravity  mode  modeled  in  section  4.4.1,  we 
observe  that  in  order  to  yield  an  averaged  (42-da.y)  rms  acoustic  travel-time  oscillation  of 
about  3  ms,  sea  surface  oscillations  of  about  5  cm  are  required.  However,  other  obser¬ 
vations  of  the  sea-surface  at  these  periods  and  record  lengths  exhibit  magnitudes  which 
are  about  ten  times  less  (Table  4.4).  Errors  from  boundary  condition  approximation  and 
discretization  cannot  account  for  this  difference.  Based  on  this  and  on  the  failure  of  the 
model  to  fit  the  data,  we  conclude  that  the  tomographic  oscillations  observed  are  not  the 
result  of  resonant  barotropic  gravity  wave  modes  of  the  Pacific  Ocean. 
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15.06  Hour  Mode 

Observed 

Model 

Receiver 

Amplitude  (ms) 

Amplitude  (ms) 

1 

2.9 

2.0 

2 

1.9 

2.2 

3 

3. 7/ 1.6 

1.6 

4 

1.3 

5 

3.3 

6 

2.4 

1.4 

7 

1.4 

1.7 

15.06  Hour  Mode 

Positive  Frequency 

Negative 

Frequency 

Receiver 

Observed  Phase 

Model  Phase 

Deviation 

Model  Phase  Deviation 

Lag  (deg.) 

Lag  (deg.) 

(deg-) 

Lag  (deg.) 

(deg.) 

1 

-161 

42 

-157 

-42 

119 

2 

-169 

29 

-162 

-29 

140 

3 

0 

0 

0 

0 

0 

4 

70 

-70 

-140 

70 

0 

5 

91/116 

-92 

177/152 

92 

1/-24 

6 

70/-86 

-97 

- 167 /- 1 1 

97 

27/- 177 

-T 

( 

-19/167 

49 

68/-144 

-49 

-30/144 

Table  4.1:  Comparison  of  acoustic  travel-time  amplitudes  (top)  and  phases  (referenced  to 
R 3 )  (bottom)  predicted  by  the  resonant  barotropic  gravity  wave  model  and  those  observed 
in  the  tomographic  data  at  the  15.06  hour  mode  (Miller  et  al.,  1988).  Phase  deviations  are 
determined  by  subtracting  the  observed  phase  lag  from  the  model  phase  lag.  Observed 
phases  are  modulated  by  ±360°  if  necessary  to  obtain  deviations  between  ±180°  (a  -175° 
deviation  is  only  10°  different  from  a  +175°  deviation).  Positive  and  negative  frequency 
indicates  clockwise  and  counterclockwise  wave  propagation.  Phase  lags  are  accurate  to 
±15°  (Appendix  13).  When  two  auto-  or  cross-spectral  peaks  occur  near  the  predicted 
mode  period,  both  are  tabulated. 
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18.46  Hour  Mode 

Observed 

Model 

Receiver 

Amplitude  (ms) 

Amplitude  (ms) 

1 

2.7 

5.4 

2 

6.1 

5,1 

3 

6.1 

6.1 

4 

5.9 

7.2 

5 

10.4 

8.7 

6 

6.5 

5.2 

7 

1.3 

4,1 

18.46  Hour  Mode 


Receiver 

Observed  Phase 
Lag  (deg.) 

Positive  Frequency 

Negative 

Frequency 

■sroHjang 

■my 

1 

-123 

8 

131 

-8 

115 

2 

150/115 

9 

-141/106 

-9 

-159/- 124 

3 

0 

0 

0 

0 

0 

4 

-120/132 

-9 

111/141 

9 

129/- 123 

5 

-6/- 154 

-15 

-9/139 

15 

21/169 

6 

126/8 

-18 

-144/-26 

18 

-108/10 

7 

69 

4 

-65 

-4 

-73 

Table  4.2:  Same  as  Table  4.1,  except  comparisons  are  made  at  the  predicted  18.46  hour 
mode. 
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21.66  Hour  Mode 

Positive  Frequency 

Negative 

Frequency 

Receiver 

Observed  Phase 

Model  Phase 

Deviation 

Model  Phase  Deviation 

Lag  (deg.) 

Lag  (deg.) 

(deg.) 

Lag  (dog.) 

(deg.) 

1 

63/174 

-4 

67/178 

4 

-59/- 170 

2 

58/71 

-2 

-60/- 73 

2 

-56/- 69 

3 

0 

0 

0 

0 

0 

4 

35/180 

5 

-30/- 175 

-5 

-40/175 

5 

53 

14 

-39 

-14 

-67 

6 

-53 

10 

63 

-10 

43 

7 

77/-107 

-8 

-85/99 

8 

-69/115 

Table  4.3:  Same  as  Table  4.1,  except  comparisons  are  made  at  the  predicted  21.66  hour 
mode. 
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FREQUENCY  (cpd) 

Figure  4-6:  Power  spectral  densities  of  42-day  time  series  of  sea  level  at  La  Jolla,  California 
(top)  and  Mokuoloe,  Oahu  (bottom).  The  La  Jolla  sea  level  data  is  from  yeardays  41  to 
83,  1981  and  the  Mokuoloe  sea  level  data  is  from  yeardays  19  to  61,  1988.  A  Bartlett 
window  is  applied  to  the  data.  The  labeled  peaks  at  periods  between  15  and  23  hours  are 
not  sidelobes  of  the  tides. 
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Fa  Jolla,  Calif,  sea-level  12-day  record 
( Figure  1-0) 

Bottom  pressure  gage  (300  km  east  of  Oahu) 
00-day  record  (Filloux,  1981) 

Bottom  pressure  gage  (1 150  km  west  of  Calif.) 
7-day  record  (Filloux,  1971) 


Fable  1.1:  A  summary  showing  amplitudes  of  sea  level  oscillations  above  background 
spectral  levels  at  periods  between  15  and  23  hours. 
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Chapter  5 


Discussion 

5.1  Summary  of  Results 

We  have  shown  in  the  previous  chapter  how  the  tomographic  data  provides  several  con¬ 
straints  on  the  problem  of  identifying  a  cause  for  the  observed  travel-time  oscillations. 
Other  experiments  in  the  Pacific  Ocean  add  further  constraints.  At  this  point,  having 
eliminated  the  more  obvious  processes  which  could  cause  these  oscillations,  it  is  useful  to 
summarize  these  constraints  and  perhaps  obtain  a  new  perspective  on  the  nature  of  an 
oceanographic  process  required  to  generate  the  observed  results.  Such  a  process  must: 

•  have  about  0.5  cm  oscillations  in  sea  level. 

•  generate,  through  current  or  temperature  perturbations,  acoustic  travel-times  of 
about  10  ms  rms  over  10-day  time  periods  and  about  3  ms  rms  over  42-day  time 
periods. 
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•  result,  in  large  plia.se  lags  in  acoustic  t  ravel-time  between  different  receivers  consistent 
with  the  data. 

•  be  intermittent  (averaging  over  time  reduces  the  magnitude  of  the  oscillations). 

•  occur  simultaneously  at  periods  between  15  and  23  hours  at  similar  magnitudes. 

The  requirement  for  small  sea  level  variations  brings  us  back  to  the  concept  of  internal 
waves.  Because  of  the  signal  processing,  our  data  is  insensitive  to  higher  order  internal 
wave  modes,  and  we  have  previously  ignored  such  processes.  There  is  a  possibility  that 
internal  waves  at  the  turning  latitude  (where  wave  frequency  equals  the  inertial  frequency) 
may  intensify  and  give  rise  to  temperature  perturbations  sufficient  to  yield  the  observed 
travel-time  magnitudes  (M.  Briscoe,  personal  communication). 

The  inertial  period  in  the  region  of  the  experiment  ranges  between  32.7  hours  (at 
Oahu)  and  about  15.7  hours  (at  R 1 ).  In  order  for  travel-time  fluctuations  to  occur  at  a 
receiver  due  to  this  process,  the  ray  path  must  cross  the  turning  latitude  for  a  particular 
period.  For  example,  at  H-5  and  R6  (inertial  period  about  24  hours),  we  expect  to  observe 
fluctuations  due  to  this  process  only  between  24  and  32.7  hour  periods.  Yet  we  observe 
prominent  oscillations  at  these  receivers  at  15  to  23  hour  periods.  Any  test  of  this  theory 
must  account  for  these  oscillations. 

For  completeness,  we  note  in  the  literature  that  several  other  theories  explaining  the 
observed  travel-time  oscillations  have  been  examined  and  rejected  (Spiesberger  et  al., 
1989b). 
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5.2  Conclusion 


YVo  must  confess  that,  vvhilo  the  results  presented  here  are  significant,  they  are  unsatisfying 
in  that  we  are  unable  to  find  an  oceanographic  process  which  explains  the  acoustic  travel- 
time  oscillations.  Had  it  not  been  for  other  work  demonstrating  agreement  between  the 
tomographic  data  and  models  of  barotropic  tides  (Headrick,  1990),  we  may  have  even 
questioned  the  validity  of  the  data  itself.  In  one  sense,  this  study  highlights  the  concept 
that  while  it  is  easy  to  design  an  experiment  to  prove  or  disprove  a  specific  theory,  it  is 
much  more  difficult  to  find  a  theory  to  fit  observed  phenomena. 

To  return  to  the  initial  motivation  for  basin-scale  acoustic  tomography,  that  of  ob¬ 
taining  synoptic  measurements  of  basin-scale  phenomena,  we  find  ourselves  both  covering 
new  ground  and  old  ideas  nomographic  measurements  provide  a  unique  window  into 
the  understanding  o'  asm-scale  processes.  Many  such  processes  have  been  theorized, 
but  have  no  direct  observational  support.  Yet  while  trying  to  understand  these  previ¬ 
ously  unobserved  phenomena,  we  can  also  obtain  a  fuller  understanding  of  “well  known” 
phenomena  (such  as  tides). 

In  this  study,  we  have  limited  ourselves  to  a  class  of  interesting  observations  in  a  very 
narrow  frequency  range.  The  data  set  however,  allows  for  observations  of  processes  having 
periods  from  about  1000  hours  to  0  minutes.  Use  of  this  data  in  conjunction  with  point 
measurements  in  the  Pacific  is  a  worthwhile  and  rich  topic  of  research. 
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Appendix  A 


Island  Trapped  Internal  Waves 


l  ho  solution  for  the  perturbed  pressure  field  of  internal  waves  scattering  from  a  cylinder 


is. 


P(r,9)  =  A  £ 


771  =  —  OO 


-  I  r^!:gyy;)  »„(*0 


r,imO 


( A.l ) 


Mirm(ka)  -  tfllm(ka ) 
where  the  primed  variables  denote  derivatives,  Jm  and  //m  are  Bessel  and  llankel  functions 
of  the  first  kind  of  order  m  respectively,  /I  is  the  amplitude  of  the  incident  plane  internal 
wave  with  frequency  a  and  horizontal  wavenumber  h\  a  is  the  cylinder  radius,  and  /  is 
the  Coriolis  frequency  (Proudman,  1914;  Wunsch,  1972).  The  origin  of  the  cylindrical 
coordinate  system,  with  radius  r  and  polar  angle  9,  is  the  center  of  the  cylinder.  The 
horizontal  wavenumber  is  related  to  the  frequency  through  the  dispersion  relationship  for 
internal  waves  by 

k1  = 


,2  _  ^  -  fl 


ghn 

where  hn  is  the  equivalent  depth  of  the  nLh  internal  wave  mode,  and  g  is  gravitational 
acceleration.  Pseudo- resonances  occur  if  proper  choices  of  k  and  a  result  in  |f*(r, 9)\  /A 
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greater  than  about  5  for  a  given  cylinder  radius. 

A  logical  choice  for  selecting  the  cylinder  radius  is  made  by  matching  the  circumference 
of  Oahu  to  that  of  a  cylinder  (29. 1  km)  (Luther,  1985).  Geographical  considerations  limit 
the  maximum  island  radius  to  10  km,  at  which  point  the  cylinder  would  touch  the 
neighboring  island  of  Molokai.  Although  not  likely,  bathymetry  suggests  the  possibility 
of  the  Oahu/Molokai/Maui  group  of  islands  being  considered  as  a  single  model  island, 
resulting  in  a  cylinder  radius  of  about  110  km. 

Historical  temperature  and  salinity  profiles  near  Oahu  are  used  to  determine  the  equiv¬ 
alent  depths  of  the  first  three  internal  wave  modes  (.88  m,  .25  m,  and  .10  m  respectively). 
We  evaluate  |  P(r,0)\  for  equivalent  depths  of  .01-1.0  m  in  .01  m  increments  and  island 
radii  of  15-50  km  in  5  km  increments  as  well  as  90-120  km  in  10  km  increments.  For  each 
possible  pairing  of  these  equivalent  depths  and  island  radii,  we  evaluate  the  magnitude 
of  the  scattering  function  between  10  and  25  hours  at  .28  hr  increments.  The  value  of 
I  /’(r,  0)|  /A  does  not  exceed  3.0  and  thus,  no  significant  pseudo-resonances  occur. 

Wo  also  search  for  poles  of  P(r,$)  in  the  complex  frequency  plane  near  the  real  fre¬ 
quency  axis  by  evaluating  a  contour  integral  (Figure  A- 1 )  for  the  same  island  radii  and 
equivalent  depth  combinations  used  in  the  previous  calculation.  These  contour  integrals 
are  all  zero,  indicating  P(r,0)  contains  no  poles  which  cause  pseudo-resonances  at  periods 
between  10  and  25  hours. 
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Figure  A- 1 :  The  contour  integral  path  of  P(r,0)  encloses  a  rectangular  area  in  the  complex 
frequency  domain,  where  the  real  axis  represents  the  domain  of  pure  oscillatory  waves. 
In  order  for  a  pole  of  P(r,0)  to  cause  a  pseudo- resonance  at  periods  between  10  and  25 
hours,  it  must  lie  within  the  region  bounded  by  Q  >  5,  where  Q  is  the  resonance  factor 
for  a  lightly  damped  oscillator  (Q  =  iRe{cr}/Im{(r}).  Points  on  the  contour  where  Q  —  5 
arc  marked  with  an  x. 
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Appendix  B 


Reliability  of  Cross- Spectral 
Phase 


We  treat  auto-  and  cross-spect  ral  peaks  as  isolated  events  which  are  too  large  to  be  caused 
by  measurement  mror  (Spiesberger  et  al,  1989b).  In  other  words,  we  are  not  interested 
in  standard  calculations  of  95%  confidence  intervals  because  we  do  not  care  if  a  certain 
acoustic  event  would  be  observed  in  95  out.  of  100  different  experiments. 

For  two  sinusoidal  signals  of  the  same  frequency  and  specified  phase  lag  (in  the  absence 
of  noise),  the  cross-spectral  phase  at  the  signal  frequency  is  equal  to  the  phase  lag  between 
the  signals.  When  dealing  with  real  data,  however,  noise  is  superimposed  causing  reduced 
coherency  and  cross- spectral  phase  errors.  Additionally,  some  form  of  averaging  must 
be  performed  on  the  cross-spectra  in  order  to  obtain  a  coherency  spectrum  (Jenkins 
and  Watts,  1968).  This  significantly  increases  the  bandwidth  of  the  cross-spectra,  with 
the  result  that  cross-spectral  phases  arc  biased  by  noise  and  other  signals  within  the 
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bandwidth  of  the  cross-sport ra. 

In  light  of  the  above,  we  estimate  the  reliability  of  cross-spectral  phases  using  the 
barotropic  tidal  signals  in  the  tomographic  data.  Headrick  (1990)  has  analyzed  the 
barotropic  tidal  signals  in  the  tomographic  data,  and  finds  close  agreement  between  the 
data  and  models  of  the  diurnal  tides.  A  continuous  J  10-day  data  segment  is  available 
from  four  receivers.  Dividing  these  data  sets  into  consecutive  55-day  sections,  we  com¬ 
pute  cross-spectra  for  each  receiver  combination  using  the  same  processing  and  averaging 
used  on  the  12-day  data  segments  discussed  in  the  body  of  this  paper.  At  the  diurnal  tide, 
K i ,  cross-spectral  phases  obtained  from  each  55-day  segment  at  each  receiver  combination 
is  consistent  within  ±15°.  We  therefore  use  ±15°  as  a  measure  of  the  reliability  of  our 
cross-spectral  phases. 
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Appendix  C 


A  Two-Mode  Gravity  Wave 
Model 


Lot  the  observed  acoustic  travel-time  at  receiver  i  be, 

c,  sin((T/  +  9t),  (C.l) 

where  c,  is  the  magnitude  of  the  rms  travel-time  estimated  from  auto-spectra,  9,  is  the 
cross-spectral  phase  (referenced  to  R3),  a  is  the  mode  frequency,  and  t  is  time.  We 
model  the  acoustic  travel-time  at  receiver  i  from  a  linear  combination  of  the  clockwise 
and  counterclockwise  travelling  resonant  barotropic  gravity  waves  (Section  4.4.1)  as, 

c,sin((7  +  9i)  =  aqflji  sin(<r  +  <ptl  +  Cl)  +  a^n^sitH-cr  +  <t>i2  +  C2)  (C.2) 

where  the  variables  with  tildes  are  the  model  estimates,  and  the  right  side  of  the  equation 
corresponds  to  the  summation  of  the  clockwise  and  counterclockwise  modes. 

Thus  for  each  mode  we  have  four  model  parameters,  (oi , oji Cl  1  k*2)*  which  are  varied 
to  minimize  the  error  between  the  observed  travel-times  and  the  model  estimates.  We 


43 


establish  a  relevant  (but  somewhat  arbitrary)  cost  function,  which  evaluates  how  well 


the  two-mode  model  fits  the  observed  data  at  all  receivers  for  a  given  mode.  After  some 
trial  and  error,  we  selected  a  cost  function  which  penalizes  percentage  errors  in  magnitude 
and  deviations  in  phase  between  the  model  and  the  data.  The  cost  function  selected  is, 

•/(Ol,OMl.<2)  =  iz  (j  +  f)  '  +  100  •  (£•  -  (C.3) 

where  the  summation  is  taken  over  all  receivers. 

For  each  mode  we  find,  after  minimizing  the  cost  function,  that  the  magnitude  and 
phase  obtained  from  the  two-mode  model  of  resonant  barolropic  gravity  waves  performed 
no  better  than  the  results  using  a  single  mode  model. 
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